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LONG PERIOD PERTURBATIONS IN THE ELEMENTS OF 
ARTIFICIAL SATELLITES DUE TO HIGH ORDER ZONAL HARMONICS 


David Fisher 


ABSTRACT 

General formulas for the long period perturbations in the elements of arti- 
ficial satellites due to high order zonal harmonics are presented. The formulas 
derived here are based on a more general determining function than Brouwer’s 
determining function in order to permit indirect gravitational effects such as 
those due to the sun and the moon to be included in the perturbations. 

The perturbations are given in the form of Fourier series in the argument 
of perigee and are associated with each cf the zonal harmonics of the earth’s 
gravitational field. The coefficients of the terms of the Fourier series are 
separable into two functions , one a power series in eccentricity and the other 
a power senes in the inclination as shown by Kaula, Reference 1. The expres- 
sion for each of these types of series is given. 

In order to simplify the calculations, recursion formulas are presented. 
These formulas show that the evaluation of the terms of the power series can be 
made in cascade fashion, each successive term being a simple multiple of the 
preceding one. 

Finally the contributions to the Lagrange planetary equations due to the long 
period even and odd zonal harmonics are given. 


iii 


RECEDING 


page blank 


*• 


n °t filmed. 


CONTENTS 


Page 


ABSTRACT iii 

INTRODUCTION 1 

THE HAMILTONIAN 4 

Secular Part of the Hamiltonian 4 

Periodic Part of Hamiltonian (n even) 5 

Periodic Part of Hamiltonian (n odd) 7 

PERTURBATIONS IN THE ELEMENTS. 9 

Long Period Perturbations - Even Zonal Harmonics 10 

Long Period Perturbations - Odd Zonal Harmonics 12 

THE MEAN MOTIONS 13 

LONG PERIOD CONTRIBUTIONS TO LAGRANGE’S PLANETARY 
EQUATIONS . 14 

Long Period Contributions To Lagrange’s Planetary Equations - 

Even Zonal Harmonics 15 

Long Period Contributions To Lagrange’s Planetary Equations - 

Odd Zonal Harmonics 16 

REFERENCES 16 

APPENDIX A-NOMENCLATURE 19 

APPENDIX B— R F CUR SION FORMULAS 21 

APPENDIX C-THE MEAN MOTION OF THE ARGUMENT OF 

PERIGEE AND SOME USEFUL DERIVATIVES 23 


v 


LONG PERIOD PERTURBATIONS IN THE ELEMENTS OF 
ARTIFICIAL SATELLITES DUE TO HIGH ORDER ZONAL HARMONICS 


INTRODUCTION 

In Reference 2, the force function for the zonal harmonics of the earth's 
gravitational field is given by 


U 


E 

r 


1 


■D 



(sin<£) 


(1) 


where /x = k 2 m e , k is the Gaussian constant, m e the mass of the earth, r the 
radius from the earth’s center to the satellite, J n the zonal coefficient, <t> the 
latitude of the satellite in equatorial geocentric coordinates. 

The definition of the Legendre Polynomials is given by 


Pn<*> 


1 d n (x 2 - l) n 
2 n n! dx n 


( 2 ) 


Brouwer, Reference 2, expanded the force function up to order 5. He then 
deduced the long period perturbations using a determining function by von Zeipel’s 
method. Kaula, Reference 1, and Giacaglia, Reference 3, extended the work of 
Brouwer to include the effects of zonal harmonics of any arbitrary order. 

It is shown in References 1 and 2 that the long period part of the Hamiltonian 
derived from the force function (1) can be represented in the form 



F(i) G(e) cos 2qg 


(3) 


1 


for even n , and 




F( i ) G(e) sin (2q + 1) g 


(4) 


for n odd. F( i) and G(e) are the inclination and eccentricity functions which 
are described below. 

In References 2 and 3 the corresponding determining functions are given by 


S 



2q 


F(i) G(e) sin 2qg 
3 M 4 J 2 b2 (* " 5 cos 2 

■7 



(5) 


for even zonal harmonics, and 



( 6 ) 


for odd zonal harmonics, where L and G are Delaunay variables with L = i\L a and 
G = L /l - e 2 , a and e are the semimajor axis and eccentricity of the satellite 
respectively. 


In this report the results of References 2 and 3 are extended by introducing 
more general determining functions than those given by Equations 5 and 6. 

Thus we define 


S 


F(i) G(e) sin 2qg 
2qa 


(7) 


2 


for even n and 


S 


ii F(i)G(e)cos(2q+l)g 


( 2q + 1 ) a 


< 8 ) 


for n odd, where 


a 



(9) 


In Equation (9) the quantity F a represente the secular part of the Hamiltonian. 
If it is desired to include indirect solar and lunar effects F a must include the 
secular part of the solar and lunar gravitational potential. The expression for 
a including solar and lunar effects as well as the zonal harmonics J 2 , J 2 a and 
J 4 is given in Appendix C. The quantity 


3 M 4 J 2 k 2 (l ” 5 cos 2 i) 

* L 3 G 4 


appearing in Equations 5 and 6 is a truncation of a and is valid for many of the 
artificial satellites which have been launched. However for satellites more 
than one or two earth radii from the surface of the earth the indirect solar and 
lunar effects may be quite significant. In such cases direct solar and lunar ef- 
fects (Reference 4) should also be included. 

In order to shorten calculations, recursion formulas are given so that simple 
multiplications are involved in generating the required series. 

The formulas given in Reference 1, are very compact and include not only 
the formulations for the zonal harmonics but formulations for tesseral harmonics 
as well. Consequently, in order to display the characteristics of the various 
components of the disturbing function, the formulations of Reference 1 are modi- 
fied and formulas are given describing separately the secular, the long period 
even, and the long period odd zonal harmonics. 
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THE 


HAMILTONIAN 


The Hamiltonian, References 2 and 3, for the high order zonals, is in this 
report divided into secular, long period even, and long period odd portions. The 
contributions to the mean motion are derived from the secular part of the 
Hamiltonian while the perturbations in the elements are derived from the long 
period parts via a determining function. 

Secular Part of the Hamiltonian 

From the results of References 2 and 3 we find that the secular part of the 
Hamiltonian is given by 




M n+ 2 J„b n 

PT 

2 n L 3 G 2n ~ 1 


( 10 ) 


where 


P = 



j a o 






T = 



B 


0,k 


sin 2k i 



To illustrate this formula we note for n « 4 


K 


o,o 


1, 



3 

1 


so that 



4 


®0,0 “ ® ’ 


B 0l | ” “ 30 # 


hence 


T - 6-30 sin 2 i + 


■y sin’ i , 


and consequently 


-^ 6 J 4 b 4 


* 4 F 16 L 3 G 


7“ ^1 + -J e 2 ^ ^6 - 30 sin 2 i + sin 4 i) 


When we write 


I k 4 • 


L 2 ’ 


• 6 • 
sin* i 


we find that this expression agrees with the result of Reference 2. 


Periodic Part of Hamiltonian (n even) 


The contribution to the periodic part of the Hamiltonian due to the n th even 
zonal harmonic J n , is given by the Fourier series 


“A* n+2 J„b n 


A F = “ 
n 2 


ix ■ u mmmrntm 

n L 3 G J"-l J2 KB cos 2qg 


where 


K = 


e 2j 
>J e 


5 



and 


B 


q.k 


o2k- I 



n + 2k\ / 2k 


2k / \k “ q , 


for the case j = 4, the Fourier series (11) consists of the single term 


-M 6 J«b 4 

^4 ^ “ 16 L 3 G 7 ^ ® cos 


where 


K 




2 


and 


B = B 1(1 sin 2 i + B l 2 sin 4 i 

B = 5 sin 2 i (6 -7 sin 2 i) 


making the substitutions for J 4 , e 2 and sin 2 k indicated in the example il- 
lustrating the secular part of the Hamiltonian we obtain agreement with Refer- 
ence 2. 
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Periodic Part of Hamiltonian (n odd) 


If we combine the results of the formulas given above, we find that the odd 
zonal harmonics define the Fourier series 


where 


A F = 

n 2 


-M"* 1 J„b n 




— Q T CDsin [(2q + 1) g] 
q*7 


( 12 ) 


■i* 


c = 


j“q 


C q,j ” 22T+T\ 


n- 1 \ /2j + 1 


2 j + 1/ \ n - q 


D = 


L 

k*q 


D ,,k sin2k+i 1 • 


and 


_ ("I) 


£^-l + q“k 


'q,k 


2 2k 


n 


n- 1 


- k 



To illustrate Equation (12), we apply it to the case n = 5. For the fifth 
zonal harmonic the Fourier series given by Equation 12 consists of the sum 


*Vj s b 5 


A S F “ 2 S L 3 G 11 S ^ n 8 + (^^q a 1 S * n * 
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We find that for q ~ 0 


C 





+ 


3 

7 e 


’) 


D - D 0 0 sin i + D 0 t sin 3 i + D 0 2 sin 5 i 


sin i 


^60-210 sin 2 i + 


315 . 

sin 


Therefore 


(CD) q= 


q=0 


e sm l 


i ^2 + e 2 j ^60- 210 sin 2 


i + 


315 

2 


sin 


■) 


For q = 1, 


C = 




D - j sin 3 i + D x 2 sin 6 i 


sin 3 



315 

~r 


sin' 


0 


So that 


(CD),,, = ( 70 . 315 s . n2 .) 


These results are in agreement with Reference 2. 
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PERTURBATIONS IN THE ELEMENTS 

The method of canonical variables is used in this report. Thus if S is a 
determining function we find from Reference 2, that 


G' 


G" + 


as 

ag' 


v 


V - 


as 

aL" 


g' 


_ _as_ 


h' 


h" 


as 

aH" 


(13) 


The determining function S is considered to be a function of the new variables 
L", G" , H" and the odd variable g' . The variables V and h' do not appear in S 
since we are dealing with long period effects in the zonal harmonics. 


We also have from Reference 2 that the mean motions of 4 , g and h are given 
by 


dV 


dt " 

dv 

dg" _ 


“ar - 

” dG" 

dh" 


dt 

aH" 


(14) 


where F s is the secular part of the Hamiltonian 


Long Period Perturbations - Even Zonal Harmonics 


We first find the determining function S by applying the definitions of Equa- 
tions (3), (4), (7) and (8) to Equations (11). Thus for n even the contribution of 
the n th harmonic to the determining function is given by 


A S n 

n w 


J n b n 
2 n L 3 G 2 "" 1 



KB sin 2qg 
2qd 


(15) 


The quantity a is the mean motion of g. In Appendix C a formula for a is 
given which includes solar and lunar effects. A more general definition of a 
for zonal harmonics is given by the second of Equation 14. 

The symbol 8 is introduced to denote the difference between an element 
and its mean value. From Equations (13) and (15) we find 


8 „G 


J„b" 
2 n L 3 G 2n ~ 1 


n- 2 



q= 1 



cos 2qg 


(16) 


Similarly we find 


8 't 
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^ n * 2 J n b n 
2 n L 3 G 2 "' 1 


n-2 



KB ( 
2qa \ e 


1 d'a\ . 

JlL ) sin 2qg 


(17) 


an expression for da/BL is given in Appendix C. The quantity is given by 


hl-L /Jl\ 

K BL 




( 18 ) 


10 



^ n * 2 Jn b " y» KB / 
2 n L 3 G 2n Z-» 2qi\ e 

a« 1 ' 



sin 2qg 


(19) 


G 2n / KB \ 
KB <?G l G 2n- 1 J 


1 - 2n 4 2c ” 8 ^kB^sin^i- ^ £ iK,.)® 21 "’ 

k*q j»q 


( 20 ) 


„ d ' a 
a o " °ac 


An expression for a is given in Appendix C. 


u n+2 J b n 
8 h = 1 Jn 

n 2 n L 3 G 2n 


n-2 


£ - i i ") sin *« 

q® i ' ' 


(21) 


where 


G dB 

' B <9H 


2 cos 




B l_, ~ “q.k 

k* q 


k B_ sin 2k ” 2 i 


( 22 ) 


a 


H 


G 


da 

<?h 5 


a H is given in Appendix C 


Long Period Perturbations - Odd Zonal Harmonics 


The determining function for the odd zonal harmonics is found by applying 
the definitions of Equations (3), (4), (7) and (8) to Equation (12). We then obtain, 
for n odd, the contribution of the n th zonal harmonic 


M n + 2 J n b n 
S ~ 2 n L 3 G 2n “ 1 


n~ 3 


CD cos [(2q + l)g] 

( 2q + 1 ) a 


(23) 


Again, by applying the formulas (13) to Equation (23) then results 


n-3 


■M n+ 2 Jn b ° A CD sin [(2q + l)g 
2 n L 3 G 2n " 1 / ^ a 


q=0 


(24) 


q* o 


where 


v. 


]J 
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±(s) 

BL \ L ») 


_1 

2 


n" 3 
~7T 


|-3 + ^l £ (2i + l)C q ,e^} 


-/‘’J.b" 


« g 2 n L 3 G 2 


q= o 


(26) 


12 




•t* 




= ~M" t2 J n b" ^ CD 
n 2 n L 3 G 2n L-i (2q + 1) 

q= 0 


I ^ 0 " \ *h) cos fc 2q + *) &] < 27 > 


G <?D 

" D <9H 


TT^ ( 2k + 1 > D q,k sin 2k ' x i 


THE MEAN MOTIONS 

The mean motions of the angular elements are given by Equations (14) 
where F f is defined by Equation (10). We then find that 






/x n+2 J n b n 

: — T“ <(1 ~ 2n) PT+ 2Pcos2 

2 n L 3 G 2n I 


n/2 


22 kB o,k sin2ii ‘ ji 


k* 1 


** 2 T(l -e 2 ) 



A 


n 



-M n * 2 J n b n 
2 "*'L 3 G 2 " 


n/2 


(cos i) p y 


k* 1 


k B 


0 ,k 


sin 


2k- 2 


(30) 


LONG PERIOD CONTRIBUTIONS TO LAGRANGE’S PLANETARY EQUATIONS 

The long period contributions of the zonal harmonics are found by forming 
the derivatives of Equations 10, 11, and 12. The derivatives of Equation (10) 
are given by Equations (28) - (30), and should be added to the corresponding 
derivatives below. It is convenient to use the following formulas for the ele- 
ments e and i 



(31) 


(32) 


The equations for the elements G , ^ , g , and h are given by the theory of 
canonical variables 


dG 

dF 

di _ 

<?F 

_ i i ' ■ « 

dt “ 

dg ’ 

dt 

dh 

dg _ 
dt 

< 9 F 

“ < 9 G ’ 

dh 

dt " 

dF 

“ aH 


( 33 ) 
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L 



Zonal Harmonics 


For n even we find 


/dG\ M n+2 J n bn ? _ , 

\\ItJ = 2L qKBsin2qg 


A 

A n Idtj 


b" p; r p; 

y - £Vc ^~ » 22 B ■ 3K+2 < i ' e,) £] jK .i 


e 2 j ' 2 1 cos 2qg 


a l dg 1 
A n Idti 


M"* 2 J n b" 
2 n L 3 G 2n ” 1 


— ^ (l - 2n) KB+ 2Xcos 2 i k B q ( 


sin 2k ' 2 i 


- 2B (l - e 2 ) j K q j e 2j “ 2 cos 2qg 


. /dh\ ^ +2 Jn b ° Tf . 
W) ~ 2 n " i" l3q2" K C ° S 1 


n-2 

T n/2 

r r 

qM k=q 


k B k sin 21 *" 2 i cos 2qg (37) 


The functions K, B, and the coefficients K q j and B q k are defined by Equa- 
tion (11). 
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Long Period Contributions 
Zonal Harmonics 



For n odd wc obtain 


/dG\ “M n+2 J n b n ^2^ r , 

A n Vdfj 2 n L 3 c j n -l 2_ ( 2C 1" 1) CD COS [( 2q + 1 ) g] 

q* 0 


n- 3 n- 3 

/dt\ . ^ 2 Jn bn ^ ^ 

'"\dtj ' 2 n L 4 G 2 "" * 2Z D_3C+ ( 1 ' e ) £ 

q“ 0 L ) * q 


(2j + l)C q j e 2i ’ 1 [sin(2q + 1) g] 


- »i~ * 

/dg\ ^ n * 2 J n b " f 3 

„(afj = 2" L 3 g~ £ yl “ 2n) CD + Ceos 2 i Y\ ( 2k + 1 ) D q ,k sinJk "' 1 

q*0 L k=q 


D (l ~ e 2 ) ^ (2j + 1) C q#i e 2 *" 1 * sin [( 2q + 1) g] 


n*3 n-1 


( au \ “ Ceos i u n+2 T b n . 2 . . 2 ■ 

*„(ar) = — 2 "l 3 g 2 " £ £ <2k + 2) D 0 ' k sin2k ‘ sin t (2q + 2) g 5 


q* 0 k* q 


The functions C and D, and the coefficients C q j and D qk are defined by Equa- 
tion (12). 
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APPENDIX A 
NOMENCLATURE 

B inclination function 

C eccentricity function 

D inclination function 

F Hamiltonian 

F a Secular part of Hamiltonian 

F( i ) inclination function 

G Delaunay variable G = l/I - e* 

G(e) eccentricity function 

H Delaunay variable H = G cos i 

J n ther» th zonal harmonic 

K eccentricity function 

L Delaunay variable L = /jla 

P eccentricity function 

P n zonal harmonic of n th degree 

S determining function 

T inclination function 

U force function 

a semimajor axis 

b mean radius of earth 

e eccentricity 

19 





g Delaunay variable - argument of perigee 
h Delaunay variable - right ascension of the node of the satellite 
i inclination of satellite orbit 

j summation index 

k summation index 

n appears as an index of the degree of the zonal harmonic 
the Delaunay mean anomaly 
q summation index 

i radius of the satellite 

a mean motion of argument of perigee 

M = k 2 m* 

k = Gaussian constant 

m ffi = mass of earth 

e c eccentricity of moon’s orbit (0.054900489) 

i ec inclination of moon's orbit plane to ecliptic ~5°.1453964 

m 0 ratio of mass of sun to combined mass of sun and earth (0.999997) 

ratio of mass of moon to combined mass of moon and earth (0.012150668) 
n Q mean motion o'' sun 0.98560027 degrees per day 
n c mean motion of moon 13.064999 degrees per day 
e obliquity of the ecliptic 
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APPENDIX B 
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When the Hamiltonian and perturbations in the elements for many zonal 
harmonics are required it is convenient to use recursion formulas to compute 
sequences of coefficients of the inclination and eccentricity functions. 

In deriving the recursion formulas the indices n, q, j and k are varied. It 
should be noted that although the coefficients in Equations 10, 11, and 12 depend 
on n this index does not appear as a subscript in order to simplify the notation. 
However when the index n is varied the additional subscript n is added to the 
coefficient as in the formulas below. 

The recursion formulas given below follow readily from the definitions in 
Equations 10, 11, and 12. 

Recursion Formulas For the Inclination and Eccentricity Functions (n even) 


(n + 1) n 

K q,j,n+2 “ (n + 1 “ 2j ) (n- 2j ) K q,i,n 


(Bl) 
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a,k,n+ 2 


4(n + l+2k) 
n + 2“2k B q.k,n 


(B2) 


- 1 (n-2j-l)(n-2j-2) 

q,j+i " 4 (j +l-q)(j +l+q) N.i 


(B3) 


-1 (n~ 2k) (n + 2k + 1) 

B q,k + i " “ 4 (k+l-q) (k+1 +q) «' k 


(B4) 


The above four formulas also apply to the secular part of the functions by 
setting q - 0. 



■■ ■’ % 


q+l. j 




J 


j + q + 1 ^q.J 


(B5) 
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Recursion Formulas For the Inclination and Eccentricity Functions (n odd 
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APPENDIX C 

THE MEAN MOTION OF THE ARGUMENT OF PERIGEE 
AND SOME USEFUL DERIVATIVES 


The motion a of the argument of perigee is given by 
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b ^ h b2 


(—S') 


3 M 4 Jj b 4 

+ IB — — 


/ G 

1 90 -192-l - 


15 ** J4 b 4 

128 q8 


128 5)^ + ( 38S + 360 


G G 2 
35 + 24 X* + 25 ^ 


G. . Gi\Hl 

L 45 L J /G 4 


21-96 ^ + (-270 + 126 ^ -~J 4 ^385- 189 ^-^ 


'8 G 2 sln £ j^ L 2 * c 2 ) |( 


"e m 0 




3/2 + (! 


n e "«* /. 1 . j • \ 

_ e ^2j3/2 \* .2 Sin l ecj 


The quantity G da/dG and denoted by a Q is 


a. 


(3 J2 b2 L 4 . 2 .\ 3 J2 2 b , 4 Li [?, 

n y — & ( 2 - 1 5 c ° s 2 w + 128 a 4 G s ^ 

+ 36 f- 25 + 48 ^ + 28 ^ cos 2 i - 30 ^154 + 132 ^- + 15 “jj 


G G 2 \ 

280- 168 l - 15° 


15 J4 b4 L 8 
128 a 4 q8 


cos 4 i 


(- 


168 + 54 


G 2 


^+36^75-28 ^ 


cos 2 i + 30 - 154 + 63 


(- 


g! 

L 2 


) 


COS 4 II 


3 ai / _3 . 2 \ 

‘8 G ^”2 sin € J 


n 2 m ' 

© © 


(l‘V) 3/ 2 
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The formula for G da/dH and denoted by i H is 





